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Abstract. The effect of antiferromagnetic (AF) superexchange (SE) interaction on dou- 
ble exchange model in two-dimensional triangular lattice is investigated in this Letter. 
The long-range ferromagnetic (FM) spin-spin correlation and strong electron-spin cor- 
relation in triangular lattice are two important factors for truncation error of FM-like 
spin-spin correlation near quarter filling by truncated polynomial expansion, which 
can be overcome by high expansion moment M. Near quarter filling, non-coplanar 
spin configuration is induced by competition between FM parallel alignment medi- 
ated by itinerant electrons and AF SE interaction between localized spins. A phase 
transition from FM to non-coplanar phase is characterized by both spin structure fac- 
tor and the change of the slope of energy in term of AF SE interaction, which is accom- 
panied by a phase separation in electronic density. At strong AF SE interaction, there is 
also another phase separation accompanying the phase transition from non-coplanar 
spin phase to 120-degree spin phase. Finally we discuss the effect of temperature on 
magnetic structure. 
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1 Introduction 

Recently noncoplanar spin configuration emerging in the ferromagnetic (FM) Kondo 
model becomes one of research hot topics in strongly correlated electrons field. It is in- 
duced by the competition between FM alignment of localized spins from Hund-interaction 
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of itinerant electron and antiferromagnetic (AF) superexchange (SE) interaction on reg- 
ular lattices and frustrated lattices [IH3|. In Ref. |1|, this kind of noncoplanar spin con- 
figuration was verified by the average of spin-spin correlation and the inner product 
of three localized spins located at one triangle in triangular lattice at quarter filling via 
classic Monte Carlo simulation combined with truncated polynomial expansion method 
(TPEM) (11151. This kind of peculiar spin configuration is stable in quite large param- 
eters space [ 3 ] . 120-degree spin configuration is easily destabilized by electron doping 
from band bottom and a noncoplanar three-sublattice ordering occurs accompanied by 
an intervening phase separation, while a noncoplanar ordering does not occur by hole 
doping from band top 0. Phase separation |6HZ) is a key factor for colossal magnetore- 
sistance effect (CMR) |SH?Q. In lightly doped region, there exists a phase separation from 
the electronic density 0. 

Geometrical frustration and /or magnetic frustration from competition between FM 
double exchange (DE) interaction and AF SE interaction arouse a broad interest, inspired 
by the discovery of peculiar transport properties, such as unconventional anomalous Hall 
effect, in some pyrochlore oxides HT0HT2H . In Ref. @, the density of states and the resis- 
tivity at various Jaf and temperature were investigated at quarter filling. In noncoplanar 
spin configuration, they are independent of the temperature. The system undergoes a 
metal-insulator transition when Jaf increases. The bad metallic behavior is consistent 
with experimental observation in frustrated itinerant magnets R2Mo207 [2J. 

In Ref. Ill, we had studied the spin-spin correlation and spin structure factor at var- 
ious electron-filling cases. However, in this Letter, we focus on the spin configuration 
around quarter filling. As shown in Ref. |T), a * quarter filling in triangular lattice, ac- 
curate FM-like spin-spin correlation can only be reproduced with higher moment M of 
polynomial expansion, e.g. M~200. Here we find that the long-range FM-like correlation 
between spin-spin correlation is responsible for truncation error of TPEM, by varying the 
dimension of systems. However, FM-like spin-spin correlation at low electron filling re- 
gion, e.g. n ~ 0.1944, can be accurately reproduced as M~ 30, as shown in Fig. 4(c) |1 1. 
This fact implies that electron-spin correlation in FM Kondo model is also important for 
truncation error of TPEM. In contrast to large truncation error of FM-like spin-spin cor- 
relation near quarter filling at low truncation moment M, most physical quantities, such 
as the energy E and the electron density n, are accurate at various parameters. 

We investigate the effect of AF SE interaction on the spin configuration near quarter 
filling. Non-coplanar spin configuration is induced by competition between FM paral- 
lel alignment mediated by itinerant electrons and AF SE interaction between localized 
spins. A phase transition from FM to non-coplanar phase is characterized by both spin 
structure factor and the change of the slope of energy in term of AF SE interaction. It is 
accompanied by a phase separation in electronic density. There is also another phase sep- 
aration accompanying the phase transition from non-coplanar spin phase to 120-degree 
spin phase at strong AF SE interaction. It was found that a phase separation also occurs 
in lightly doped case in Ref. [3|. Finally we discuss the effect of temperature on magnetic 
structure. 



3 



2 Hamiltonian, triangular lattice and spin structure factor 

The hamiltonian consists of two parts, i.e., H = Hde+Haf- Double-exchange model Hoe 
describes the Hund's-rule interaction between itinerant electron and localized spins, 

H DE = -t £ {Cl cl C jA +hc)-] H Y J Cl K a afi C iifi -S i . (2.1) 

<ij>A i,a,j8 

Here the nearest-neighbor hopping integral t is adopted as the energy unit and Jh is the 
Hund interaction strength. C\ a (Cj A ) creates (annihilates) one electron at site i with spin 
dl, < ij > stands for the nearest neighbors of lattice site, and cr a a is the Pauli matrix. The 
localized spin S; at site i is assumed as 1 here. H^f = jAF*L<ij>Si-Sj describes AF SE 
interaction between nearest-neighboring lattice sites. Localized spins is treated as a clas- 
sic field (p and electron degree of freedom can be integrated for any given localized spin 
configuration. Classic Monte Carlo simulation combined with TPEM |4,5| is adopted 
to simulate the energy, electron density, spin-spin correlation and spin structure factor. 
There are three key parameters for TPEM, i.e., M, e p and which control the accu- 
racy and computational speed (Tl|4l|5|. The maximal Monte Carlo step is 40000, and the 
physical quantity is evaluated every 20 steps after first 6000 warmup steps. 

6x6 triangular lattice is shown in Fig. [T{a) and the lattice constant is denoted as a. 
Periodic boundary condition is adopted in simulation. Each lattice site has six nearest 

neighboring lattice sites. Two basic lattice vectors are U\ = (1,0) a and a^ = (%,^k)a, as 
shown in Fig. []Jb). The reciprocal lattices are b\ = -^{^, — \) and £>2 = ^^(0/1)/ an d 
the first Brillouin zone is also shown in Fig. [He). For LxL triangular lattice, similar to 
Ref. 1 1 1, the momentum q = j-bi + ^$2 is shortened as ((71,172)/ with 171 = ™, 172 = \, and m, 
n as integers from to L. Spin structure factor S(q) is defined by 

S(q) = -^E< S rSj>e i ^, (2.2) 

where < S, • Sj > and r,y is the average spin-spin correlation and the displacement between 
site i and j respectively, and N = LxL. Here j = i x + (L — 1) * L is the ordering of lattice site 
with indices (i x Jy), in which 1 < i x < L and 1 <L < L. Due to periodic boundary condition, 
only < Si • Sj > is shown. 

3 Result and discussion 

3.1 The truncation error of TPEM related to the range of FM spin-spin corre- 
lation and electron-spin correlation 

It is well known that the truncation error decreases when the truncation moment M of 
TPEM increases. Therefore the physical quantity can be exactly reproduced when M 
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approaches infinity. As discussed in Ref. [I], a large truncation moment (M > 200) is 
required to reproduce FM-type spin-spin correlation < Si • S; > in 2D triangular lattice at 
quarter-filling case, as shown in Fig. |2]T). Since there is no stable FM phase for double 
exchange model in these systems with periodic boundary condition at quarter filling, 
FM spin-spin correlation is short-ranged. Therefore, M ~ 30 and M ~ 50 is enough to 
reproduce the spin-spin correlation < Si • S; > at quarter filling in ID and 2D regular lattice 
respectively as shown in Fig. |2ja,c). At strong AF SE interaction, e.g., Jaf = 0.1, spin-spin 
correlation is FM and even AF in short range, M = 30 is enough as shown in Fig. |2fb,d) 
and Fig. 4(e) in Ref. (TJ. For 2D square lattice, a spin-flux phase II13II14I occurs, where four 
localized spins within each square lie (anti)clockwise at the same plane. For 2D triangular 
lattice, a new type of spin-flux phase, i.e., noncoplanar spin configuration, occurs as a 
result of competition of electron-mediated FM alignment and AF SE interaction |l]-[3]. 
These results show that the range of FM-type spin-spin correlation and thus FM magnetic 
domain is an important factor for the accuracy of spin-spin correlation by TPEM. 

On the other hand, we explore the effect of electron mediation on the accuracy of spin- 
spin correlation by TPEM. In three-dimensional (3D) cube lattices, the magnetization at 
quarter filling is accurate at M = 30 (3, compared with exact diagonalization method. 
Though FM-type spin-spin correlation at low electron filling region (e.g., n=0.1944) is 
long-ranged in 2D triangular lattice, FM-type spin-spin correlation is accurate at M = 
30, as shown in Fig. 4(c) (H. The electron spin correlations in both cases in above are 
weaker than that in triangular lattice near quarter filling, therefore electron mediation 
and electron-spin correlation is the other important factor for the accuracy of spin-spin 
correlation by TPEM. 

Compared with spin-spin correlation, the energy can be reproduced by TPEM at 
M~50, irrelevant of the strength of AF SE interaction Jaf , as shown in Fig. [2](f ) . Therefore 
most of physical properties are correct by TPEM, even if long-range FM spin-spin corre- 
lation in FM phase is overestimated. We notice that the energy first increases, reaches 
a maximum and then decreases, as Jaf increases from to 0.1, which indicates a phase 
transition from FM into non-coplanar spin configuration, to be discussed in the follow- 
ing. 

3.2 Effect of AF superexchange interaction on spin configuration 

Due to simplicity, one triangular lattice is chosen to show the spin configuration under 
various AF superexchange interaction. The possible electron density n e of one triangular 
lattice is 0, 1/3, 2/3 and 1. The spins are parallel at n e = 1/3 at weak AF SE interac- 
tion, while they are in 120° at n e = 2/3 and 1, independent of the value of Jaf ■ When 
Jaf is 0.17~0.18, the electron density sharply changes from 1/3 to 2/3, which indicates 
a phase separation, as shown in Fig. |3ja). The energy E first increases, reaches a max- 
imum -0.49071 and then monotonically decreases, when Jaf increases from to 0.4, as 
shown in Fig. |3lb). Interestingly, the maximal energy E max occurs in the vicinity of the 
phase separation. Meanwhile the spin-spin correlation changes from 1 to -0.5, when Jaf 
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increases from to 0.4, as shown in Fig. |3{c). The spin configuration evolves from parallel 
to non-coplanar configuration to 120°. The probability distribution of the inner product 
of the nearest-neighboring localized spins S\-S2 is analyzed at Jaf = 0.1, 0.17 and 0.4, as 
shown in Fig. Eld). The mean value of (Si • S 2 ) is 0.8381, 0.1753 and -0.48842 respectively. 
It clearly shows that the peak of probability distribution is located at 1, -0.5 for FM phase 
and 120° . However, the probability distribution has two peaks at 1 and -0.5 respectively 
for non-coplanar spin configuration. We stress that the phase separation in one triangular 
lattice is different from those in large triangular lattice, since n e = 2/3 induces 120° spin 
configuration. Therefore non-coplanar spin configuration in one triangular lattice exists 
in a small space of parameters. 

Here we investigate the effect of AF superexchange interaction on the energy E, the 
electron density n, the spin-spin correlation (Si-Sy) and spin structure factor S(q), by 
continuously tuning the strength of AF superexchange interaction Jaf- As Jaf increases 
from to 0.038f, the energy E increases, as shown in Fig. H|a). At small Jaf, the spin con- 
figuration is in FM phase and (S, • Sy) is positive. The FM phase is strengthened by the six 
nearest neighboring lattice sites in 2D triangular lattice [1 j, compared with unstable FM 
phase in 2D square lattice. When Jaf increases, FM-like spin-spin correlation decreases 
and the angle between nearest neighboring localized spins approaches to 90°, as shown 
in Fig. Htb). As a result, the spin structure factor at the vector (0,0), i.e. S(0,0), decreases 
as Jaf increases, as shown in Fig. SJc). Since E<,y>(S ! -Sy) remains positive, the energy E 
increases with Jaf by a competition between increasing Jaf and decreasing (S, • Sy) . 

As Jaf is larger than a critical value J c A l F = 0.038f, the energy E decreases from a max- 
imum. Therefore the slope of the energy in term of Jaf changes sign from positive to 
negative, as shown in Fig. Ufa). The change in the slope of the energy usually indicates 
a phase transition |T5|. Interestingly, a phase separation occurs around this phase tran- 
sition, which is indicated by the electron density, as shown in Fig. |Ua). As Jaf > Iaf> 
the energy corresponding to FM phase is too large and a new phase with a lower en- 
ergy emerges, i.e., noncoplanar spin configuration. In this phase, the spin-spin correla- 
tion between any two nearest neighboring localized spins is negative, as shown in Fig. 
Hpb), while (Si x Sy) -S^ with the sites i, j and k belonging to one triangle is nonzero |T). 
Noncoplanar spin configuration is different from both FM phase and 120° configuration 
phase. The spin structure factor S(q) has three peaks, which is located at vectors (0.5,0), 
(0,0.5) and (0.5,0.5) respectively, as shown in Fig. |4{c), which is consistent with Ref. H2J. 
Rescaled spin structure factor S(q) / N~0.29 is close to its maximum value 1/3, and finite 
size effect is small HTJ. 

When Jaf > Jaf = 0.225f, the system eventually evolves into 120° spin configuration 
phase. In this phase, the spin-spin correlation between any two nearest neighboring lo- 
calized spins closes to —0.5, as shown in corrected Fig.4(a,d) [H, while (S, x Sy) • with 
the sites i, j and k belonging to one triangle is zero (TJ. The spin structure factor S(q) has 
two peaks, which is located at vectors (1/3,2/3) and (2/3,1/3) respectively, as shown 
in corrected Fig. 5(a,d) dED, which is consistent with Ref. \2\. Rescaled spin structure 
factor S(q)/N is close to its maximum value 1/2 Q]|2). There is another phase separation 
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at J C A 2 F = 0.225f indicated by the electron density. Interestingly this phase separation in- 
dicates phase transition from noncoplanar spin configuration to 120-degree spin phase, 
which is consistent with the phase diagram shown in Ref. [2| in the limit of infinite Hund 
interaction. 

3.3 Effect of temperature on spin configuration 

We further investigate the effect of temperature on the spin-spin correlation at Jaf = 
and 0.1 respectively as shown in Fig. |5ta,b). As the temperature T increases from to 
O.lf, FM-type spin-spin correlation at Jaf = decreases as shown in Fig. |5{a), and the 
spin configuration tends to be paramagnetic at high temperature. On the other hand, 
the absolute value of the spin-spin correlation in non-coplanar spin configuration also 
decreases as T increases, as shown in Fig. [Sjb). When Jaf increases from to 0.1, the peak 
of spin structure factor moves from vector (0,0) to vectors (0,0.5), (0.5,0) and (0.5,0.5). 
Meanwhile S(0,0) decreases and S (0,0.5) increases. The rescaled peak of spin structure 
factor S(q x ,q y ) decreases when the temperature T increases, as shown in Fig. |5lc). All 
rescaled peaks collapse at T = 0.1, which implies the system is in paramagnetic phase. 
The dependence of spin-spin correlation on temperature and no phase transition at high 
temperature is consistent with the phase diagram shown in Ref. (2|. 

4 Conclusion 

We investigate the effect of AF SE interaction on the spin configuration near quarter fill- 
ing and discuss the origin of the truncation error of long-range FM-like spin-spin correla- 
tion by truncated polynomial expansion. Non-coplanar spin configuration is induced by 
competition between FM parallel alignment mediated by itinerant electrons and AF SE 
interaction between localized spins. A phase transition from FM to non-coplanar phase is 
characterized by both spin structure factor and the change of the slope of energy in term 
of AF SE. It is accompanied by a phase separation in electronic density. At strong AF 
SE interaction, there is also another phase separation accompanying the phase transition 
from non-coplanar spin phase to 120-degree spin phase. Finally we discuss the effect of 
temperature on magnetic structure. 
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Figure 1: 6x6 triangular lattice (a), lattice vector (b), and reciprocal lattice and the first Brillouin zone (c). 
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Figure 2: The comparison of spin-spin correlation <S\-Sj> between exact diagonalization (ED) and TPEM 
in (a,b) one-dimensional chain with L = 64, (c,d) two-dimensional 6x6 square lattices and (e) 6x6 triangular 
lattice. Jap = in (a,c,e) and J^p = 0.1 in (b,d). The case for J^p = 0.1 in 6x6 triangular lattice had been 
shown in Fig. 4(e) [I], (f) The relation of the energy £ to the strength of antiferromagnetic superexchange 
interaction J AF calculated by ED and TPEM with M = 30 and M = 50. T = 0.03f, fe = 8, ]4 = -8, e p = 10~ 5 
and <?,,. = 10~ 3 . (e) e p = \Q- 12 . 
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Figure 3: The electron density n (a), the energy E (b), the spin-spin correlation <S;-S;> (c), and the probability 
distribution of S\-S2 (d) at various ]af- The structure is one triangular lattice. Other parameters are T=0.01f, 
7h = 8 and }i = — 8. 




Figure 4: The energy E and the electron density n (a), the spin-spin correlation <Sj Sy> (b), and the spin 
structure factor S{q x ,cjy) (c) when } A p changes from to 0.1. T = 0.03f, Jh = 8 and ^ = —8. 
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Figure 5: The dependence of spin-spin correlation (a) <S\-Sj> at /^f = and (b) /^f =0.1 on the temperature, 
and (c) the dependence of rescaled peak of S{q x ,qy) on the temperature. 



